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Abstract — This paper considers the use of energy harvesters, 
instead of conventional time-invariant energy sources, in wireless 
cooperative communication. For the purpose of exposition, we 
study the classic three-node Gaussian relay channel with decode- 
and-forward (DF) relaying, in which the source and relay nodes 
transmit with power drawn from energy-harvesting (EH) sources. 
Assuming a deterministic EH model under which the energy 
arrival time and the harvested amount are known prior to 
transmission, the throughput maximization problem over a finite 
horizon of N transmission blocks is investigated. In particular, 
two types of data traffic with different delay constraints are 
considered: delay-constrained (DC) traffic (for which only one- 
block decoding delay is allowed at the destination) and no-delay- 
constrained (NDC) traffic (for which arbitrary decoding delay 
up to N blocks is allowed). For the DC case, we show that the 
joint source and relay power allocation over time is necessary 
to achieve the maximum throughput, and propose an efficient 
algorithm to compute the optimal power profiles. For the NDC 
case, although the throughput maximization problem is non- 
convex, we prove the optimality of a separation principle for 
the source and relay power allocation problems, based upon 
which a two-stage power allocation algorithm is developed to 
obtain the optimal source and relay power profiles separately. 
Furthermore, we compare the DC and NDC cases, and obtain 
the sufficient and necessary conditions under which the NDC case 
performs strictly better than the DC case. It is shown that NDC 
transmission is able to exploit a new form of diversity arising 
from the independent source and relay energy availability over 
time in cooperative communication, termed "energy diversity", 
even with time-invariant channels. 

Index Terms — Energy harvesting, relay channel, decode and 
forward (DF), cooperative communication, energy diversity. 

I. Introduction 

IN conventional energy-constrained wireless communication 
systems such as wireless sensor networks (WSNs), sensors 
are equipped with fixed energy supply devices, e.g., batteries, 
which have limited operation time. When thousands of sensors 
are deployed in a hostile or toxic environment, recharging or 
replacing batteries becomes inconvenient and even impossible. 
Hence, harvesting energy from the environment is a much 
easier and safer way to provide almost unlimited energy 
supply for WSNs. However, compared with conventional time- 
invariant energy sources, energy replenished by harvesters 
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is intermittent over time, e.g., energy fluctuation caused by 
time-dependent solar and wind patterns. As a result, wireless 
devices powered by renewable energy are subject to the 
energy-harvesting (EH) constraints over time, i.e., the total 
energy consumed up to any time must be less than the energy 
harvested by that time. 

Wireless communication with EH nodes has recently drawn 
significant research attention. In (TJ, 13, the authors investi- 
gated the power management strategies for WSNs with EH 
nodes, for which random EH models were assumed. For the 
point-to-point communication powered by EH sources, the 
power management problem was studied in Q, |4[ with the 
deterministic EH model, and in (5) — (Vj with the random EH 
model. In particular, with the deterministic EH model, under 
which the energy amount and arrival time are assumed to be 
known prior to transmission, the authors in Q studied the 
throughput maximization and transmission time minimization 
problems for the point-to-point additive white Gaussian noise 
(AWGN) channel. These results were generalized in Jffl by 
further considering the finite energy storage limit. With the 
block Markov random EH model, the authors in studied 
the throughput maximization problem over the fading AWGN 
channel and derived the optimal power allocation polices via 
dynamic programming and convex optimization techniques. In 
addition, with an independent and identically distributed (i.i.d.) 
EH model, the authors in [7;] studied the AWGN channel 
capacity under the EH constraints, and showed that even with 
the time-varying energy source, the same capacity can be 
achieved as that for the conventional case of constant power 
supply with the same total transmission energy consumed. It is 
worth noting that the authors in JU considered the throughput 
maximization problem for the Gaussian two-hop relay channel 
without considering the direct link between the source and the 
destination, which is a special case for the relay channel model 
considered in our paper. 

On the other hand, node cooperation has been known as an 
effective way to improve the system capacity and diversity 
performance in wireless networks. With conventional time- 
invariant energy sources, the full-duplex relay channel has 
been thoroughly investigated in, e.g., I9l- fl2l . where various 
achievable rates with decode-and-forward (DF) and compress- 
and-forward (CF) relaying schemes were obtained. For the 
half-duplex relay channel in which the relay needs to transmit 
and receive over orthogonal time slots or frequency bands, 
the achievable rates and power allocation polices have been 
examined in lfl3ll . In particular, the orthogonal half-duplex 
relay channel, in which the relay-destination link is orthogonal 
to the source-relay and source-destination links, has been 
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studied in fBl , 

In this paper, we study the half-duplex orthogonal Gaussian 
relay channel with EH source and relay nodes, as shown in Fig. 
Q] It is assumed that the relay transmits and receives over two 
different frequency bands, and thus the relay-destination link 
is orthogonal to both the source-relay and source-destination 
links. Here we consider the simple case with deterministic 
source and relay energy profiles, corresponding to practical 
scenarios where the EH level can be predicted with negligible 
errors, and leave the more general random cases for future 
study. Moreover, we focus on the DF relaying scheme for the 
purpose of exposition. We examine the throughput maximiza- 
tion problems over a finite horizon of TV-block transmission^ 
In each block, the source transmits a new message, which 
is received and decoded by the relay, and then forwarded to 
the destination in the subsequent one or more blocks. Our 
main objective is to study the structure of the optimal power 
and rate allocation at the source and the relay over different 
blocks to maximize the total throughput, under individual 
source and relay EH constraints. Specifically, we consider the 
following two types of data traffic with different decoding 
delay requirements at the destination: 

1) Delay-constrained (DC) traffic: The destination is re- 
quired to decode the i-th source message, i = 1, • • ■ , TV, 
immediately after it receives the signals from the source 
in the i-th block and from the relay in the (i+l)-th block. 
With such a requirement, the relay needs to forward the 
source message received in one block to the destination 
immediately in the next block; 

2) No-delay-constrained (NDC) traffic: The destination can 
tolerate arbitrary decoding delays provided that all source 
messages are decoded at the end of each TV-block trans- 
mission. Consequently, the relay is allowed to store the 
decoded source message of the i-th block, and forward it 
to the destination in any of the remaining (i + l)-th, ■ • • , 
(TV + l)-th blocks. 

Clearly, the NDC case allows more flexible relay operations 
than the DC case, and is thus expected to achieve a larger 
throughput in general. It is worth noting that in practical 
EH systems, the source and relay may have independent 
energy arrivals over time; as a result, there exists a new form 
of diversity, termed "energy diversity", to be exploited in 
cooperative communication with EH nodes, even for the case 
of time-invariant channels. As will be shown in this paper, the 
NDC transmission is able to exploit the energy diversity to 
achieve a larger throughput than the DC counterpart, thanks 
to the more relaxed delay requirement. The main contributions 
of this paper are summarized as follows: 
1) For the DC case, we formulate a convex throughput 
maximization problem, and develop a joint source and re- 
lay optimal power allocation algorithm by exploiting the 
Karush-Kuhn-Tucker (KKT) optimality conditions and 
the monotonic property of the optimal power allocation 
that is non-decreasing over time. It is shown that the 
developed algorithm is a forward search in time over 

'Note that in total (N + l)-block time is needed for each ./V-block 
transmission due to the one-block decoding delay at the relay. 
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Fig. 1. Orthogonal relay channel with energy harvesting source and relay 
nodes. 

the two-dimensional (at both source and relay) harvested 
energy profiles, which can be considered as an extension 
of the one-dimensional (at source only) search algorithm 
in (3), 121 for the case of point-to-point AWGN channel. 
2) For the NDC case, although the throughput maximization 
problem is in general non-convex, a separation principle 
for the source and relay power allocation problem is 
proved to be optimal, upon which the original problem is 
decoupled into two convex subproblems that separately 
achieve the optimal source and relay power allocation. 
Such optimal source and relay power allocation is shown 
to be also non-decreasing over time, similar to the DC 
case. Moreover, we derive the necessary and sufficient 
conditions for the NDC case to strictly improve the 
system throughput over the DC case. 
The remainder of the paper is organized as follows. Section 
II presents the system model and summarizes the main as- 
sumptions in this paper. Section III formulates the throughput 
maximization problems for the DC and NDC cases, respec- 
tively. Sections IV and V develop algorithms to solve the 
formulated problems for the DC and NDC cases, respectively. 
Numerical results are presented in Section VI to validate the 
theoretical results. Finally, Section VII concludes the paper. 

Notation: log(-) and ln(-) stand for the base-2 and natural 
logarithms, respectively; C(x) = | log (1 + x) denotes for the 
AWGN channel capacity; min {x, y} and max {x, y} denote 
the minimum and maximum between two real numbers x and 
y, respectively; (x) + = max(0,a;). 

II. System Model 

We consider the classic three-node relay channel, which 
consists of one source-destination pair and one relay, as shown 
in Fig. Q] We assume that the relay node operates in a half- 
duplex mode over two orthogonal frequency bands, while 
the source-relay and source-destination use the same band. 
For simplicity, we do not consider the bandwidth allocation 
problem for the relay, and assume that the source-relay and 
relay-destination links operate with equal bandwidth. 

We consider the DF relaying scheme, which requires the 
relay to successfully decode the source message. Moreover, we 
adopt an TV-block transmission protocol described as follows: 
During each of the TV source transmission blocks, say, the i- 
th block, 1 < i < TV, the source transmits a new message 
Wi with power Pg(i) and rate i?(i); upon receiving the signal 
from the source, the relay decodes Wi, and generates a binning 
index for based on the "random binning" technique |[l4l 
with rate i?s(i + l). In the (i + l)-th block, the relay transmits 
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a message Uj+i with power P R (i + 1) and rate C {Pr{i + 1)). 
It is noted that for the DC case (defined in Section I), Vi+i 
is the binning index of source message iUj only; while for 
the NDC case (defined in Section I), Vi+i may contain the 
information of binning indices for all source messages Wk's, 
k < i. Moreover, we assume that each block has B channel 
uses, where B is assumed large enough such that the channel 
capacity results in (9), lfl4l are good approximations to the 
communication rates in practical systems. 

In addition to the block transmission model, we assume that 
the harvested energy arrives at the beginning of each block 
with known amounts Es(i) in the i-th block and Er(i + 1) 
in the (i + l)-th block, i = 1, 2, • • • , N, at the source and the 
relay, respectively. In this paper, we assume that the battery 
capacity to store the harvested energy is infinite, and the 
consumed energy at the source or relay other than transmission 
energy is small and thus negligible. Thus, the amount of energy 
available for each block transmission is constrained by the 
following source and relay EH constraints: 

Z) p 5(i)<4X)^(i). k = l,— ,N, (1) 

?:=i i=i 

^Pii(i + 1) < ^^^(i + l), k = l,---,N. (2) 

i=l i=l 

For the i-th source and the (i + l)-th relay transmission 
blockH i = 1, • • • , JV, the channel input-output relationships 
are given as: 

y sr (i) = \/JhrX s (i) + n r (i), (3) 
Vsd{i) = \fh~ d x s {i) + n d (i), (4) 
Vrd{i + 1) = \J~K~dX r {i + 1) + w d {i + 1), (5) 

where x s (i) and x r (i + 1) are the transmitted signals in the 
i-th source and the (i + l)-th relay transmission blocks with 
power Ps(i) and Pn(i+1), respectively; y sr (i) is the received 
signal at the relay; y sd (i) and y r d("i+^) are the received signals 
at the destination from the source and the relay, respectively; 
h sr , h r d, and h s d are the constant channel power gains for the 
source-relay, relay-destination, and source-destination links, 
respectively; n r (i), rid(i), and Wd(i + 1) are i.i.d. circularly 
symmetric complex Gaussian (CSCG) noises each with zero 
mean and unit variance. 

With the above assumptions, the received signal-to-noise 
ratio (SNR) for the source-relay, source-destination, and relay- 
destination links are given as 7 sr (i) = Ps(i)h sr , "f s d(i) = 
Ps{i)h s d, and ^ rd (i + 1) = P R (i + l)h rd , respectively. 

2 Note that the (i + l)-th relay transmission block in fact corresponds to 
the i-th source message in the DC case. 



Define new source/relay energy and power profiles as Es(i) = 
E s (i)h sr , E R (i + 1) = E R (i + l)h rd , P s (i) = P S (i)h sr , 
and Pr{i + 1) = Pr(i + l)h r d, and new channel gains as 
h sr = h r d = 1 and h s d = jf^- = ^-o- It is easy to check that 
with the new defined parameters, each link has the same SNR 
as before and the EH constraints given in ([T]i-(|2]i are satisfied 
accordingly with the new power and energy profiles. As such, 
we could always determine the source and relay power profiles 
Ps(i)'s and Pr{i + l)'s first under the assumptions that 
h sr = h r d = 1 and h sd = jf^, and then obtain Pg(i)'s and 
PR(i + l)'s by scaling accordingly. Therefore, for notation 
simplicity and without loss of generality, we simplify the 
signal models in ©-(HJ as 

y sr (i) = x s (i) + n r (i), (6) 

Vad(i) = Vh)X s (i) + n d {i), (7) 

!/rd(* + l) =x r {i + l)+w d {i + l), (8) 

by setting h sr = h rd = 1 and h s d = ho. 

Moreover, it is assumed that < ho < 1, which means that 
the source-relay link is stronger than the source-destination 
link. Thus, the relay can always help with increasing the 
achievable rate from the source to the destination. 

III. Problem Formulation 

A. Delay-Constrained Case 

First consider the DC case. Since in the i-th source trans- 
mission block, the source transmits message Wi with power 
Ps{i) and rate R(i), the relay decodes Wi reliably only if 

R(i)<C(P s (i)). (9) 

Then in the next block, the relay partitions w; into bins with 
an equivalent rate + 1) lfT4l . and transmits the binning 

index in message i>j + i to the destination with power Pr(i + 1). 
At the destination, it first decodes if 

R B (i + l)<C(P R (i + l)), (10) 

and then decodes the original message Wi if 

R(i) <C{h P s (i)) + R B {t + l) 

<C{h P s (i)) + C(P R (i + l)), (11) 

where the second inequality is due to ( TTOb . From (O and dTTb , 
the achievable rate for the i-th source message is given by 

= min {C (P s (i)),C(h P s (i)) + C (P R (i + 1))} . 

(12) 

Note that for the case of h = 0, the coding scheme for the 
relay can be simplified to repetition coding, i.e., the source 
and the relay can use the same codebook. 



1 - 

(PI) max 9 ,„ x1 T mh{f: (P s (i)) ,C (h P s (i)) + C (P R (i + 1))} (13) 



1=1 



s. t. ©, ©, P s (i) > 0, P R (i + 1) > 0, i = 1, • ■ • , N. (14) 
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Considering the TV -block transmission, the average through- 
put in the unit of bits/sec/Hz (bps/Hz) is maximized by solving 
Problem (PI) in (fl~3T> - (TT4l> , where the factor \ in ( fT3T > is due 
to half-duplex relaying, and -jA^ is due to the fact that each 
TV -block transmission requires (N + l)-block duration. Next, 
some properties of the optimal power allocation solution for 
Problem (PI) are revealed. 

Proposition 3.1: For < ho < 1, there exist optimal power 
profiles Pg(i)'s and P R (i + l)'s, which achieve the maximum 
throughput of Problem (PI) and satisfy the following inequal- 
ities: 

C (P* s (i)) > C (WKO) + C (P* R (i + 1)) , i = 1, ■ • • ,N; 

(15) 

moreover, for the case of ho = 0, there exist optimal power 
profiles for Problem (PI) with 

= + *=!,••• ,N. (16) 

Proof: If ( fT3T > is not satisfied for any i, we can always 
decrease P R (i + 1) until it is satisfied, without reducing the 
achievable rate of the z-th source message. Similarly, if ( fT6] l is 
not satisfied at any i, we can reduce PJ(£) if -PJ(i) > ^r(*+1) 
or P R (i + 1) if P R (i + 1) > Pg(i) until the equality holds, 
without reducing the rate of the i-th source message. Thus, 
this proposition is proved. ■ 

Remark 3.1: From Proposition 13.11 we infer that the opti- 
mal solution of Problem (PI) is not unique in general, e.g., 
when the energy harvested at the relay is excessively larger 
than that at the source. In the sequel, we are only interested 
in finding the optimal solutions for Problem (PI) satisfying 
< fT~5T > and ( fT6b for the cases of < ho < 1 and ho = 0, 
respectively, which achieve the minimum energy consumptions 
at the source and relay. 

By ( TBI ) and ( fTUb . we obtain that 

C(P* s {i))>C(hoP* s (i))+C(P* R (i + l)) 

>C(h P* s (i)) + R* B (i + l), i=l,---,N. 

Together with ( [Tol l, it follows that 

R B (i + 1) = min {C{P R {i + l)),C(P s (i)) - C(h P s (i))} , 

i = l,---,N. (17) 

As such, if we can solve Problem (PI), by further applying 
(fTTT i, we can obtain the optimal relay rate allocation for the 
DC case. 

It is easy to verify that Problem (PI) is convex, and thus 
solvable by existing convex optimization techniques, e.g., the 
interior point method ifTBI . However, such an approach does 
not reveal any insight for the optimal solution. Thus, in this 



paper, we develop an alternative method to solve Problem (PI) 
by exploiting its special structure, as will be shown later in 
Section IV. 

B. No-Delay-Constrained Case 

For the NDC case, the relay operates the same as the DC 
case except that it is allowed to transmit the binning index for 
message Wi in messages Uj+i, ■ • • , ujv+i instead of Vi + i only 
as in the DC case. At the destination, the binning indices for 
all source messages can be successfully decoded if 

N N 

5^i2 fl (i + l) = 5^C(i^(i + l)), (18) 

i=l i=l 
N N 

Y,Rb(i + 1) <Y,C(P R (i + l)), 2<k<N, (19) 

i—k i—k 

which is equivalent to 

fe k 

+ >^C(P fl (i + l)), 1 < fc < JV - 1, 

i=l i=l 
N N 

5^iiB(i+l)=53C(Pfl(i+l)). (20) 

i=l i=l 

With the decoded binning index, the i-th source message 
can be decoded successfully at the destination if R(i) < 
C{h Ps(i)) + Ps(« + 1), i = l, ■•■ ,N. Combining this with 
(O, the achievable rate of the i-th source message is given as 

R(i) = mm{C(Ps{t))X(h Ps{i)) + RB(i + l)} 

= C(h P s (i)) + R B (i + 1), i = l,---,N, (21) 

where the second equality is due to a similar argument as 
Proposition 13. II Note that for the special case of ho = 0, we 
have R(i) = Rs{i + 1) in d2Tb . In addition, d2TT i implies that 

C{h P s (i)) + R B (i + 1) < C{P s {i)), i = l,--- ,N, which 
leads to 

fe k 

^C(PsW) - C(h P s (i)) >Y, R s(i + 1), k = l,---,N. 

i=l i=l 

(22) 

From (f20b and ( l22b . we obtain 

fe 

J2C(hoPs(i)) + C(P R (i + l)) 

i=l 

k 

<Y, c ( p s(i)), k = l,--- ,N. (23) 

i=l 



(P2) { P sW} X(, + m WTT) g ° {hoPs{l)) + C {Pr{i + 1)} (24) 

fe 

s. L^2C{hoP 8 (i)) + C (P R (i + 1)) < J2 C ( p s(i)) , k = 1, • • • , and (25) 
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Using ( |2"TT i and j23l , the average throughput for the NDC 
case is maximized by solving Problem (P2), which is non- 
convex due to the first constraint in (|25T >. and thus difficult 
to solve at a first glance. We will derive the optimal solution 
for this problem based on a separate source and relay power 
allocation strategy in Section V. We conclude this section by 
the following two propositions regarding Problems (PI) and 
(P2). 

Proposition 3.2: The maximum value of Problem (P2) is 
no smaller than that of Problem (PI). 

Proof: By Proposition 13.11 it follows that Problem (PI) 
has the same maximum value as Problem (PI*) in (f26b-(f27T). 
It is easy to see that any solution that satisfies d27l ) will also 
satisfy d25l l for Problem (P2), but not vice versa. As such, the 
feasible set of Problem (P2) contains that of Problem (PI*), 
which implies that the maximum value of Problem (P2) is no 
smaller than that of Problem (PI*) or (PI). The proof is thus 
completed. ■ 

Proposition 3.3: For any optimal power profiles Pg(i)'s 
and P R (i + l)'s for Problem (PI) or (P2), if < ho < 1, 
the constraint P^iLi -Pg(i) < 2~2iLi Es(i) must be satisfied 
with equality; if ho = 0, at least one of two constraints 
BEtiPSd) < EtiEsii) and BEf =1 P^+l) < 
2~2i=i pR,{i + 1) must be satisfied with equality. 

Proof: Supposing that the above equalities do not hold, 
we can increase the source or relay transmission power without 
violating the energy constraints to further improve the through- 
put, which contradicts the fact that Pg(i)'s and P R (i + l)'s 
are optimal solutions of Problem (PI) or (P2). Thus, this 
proposition is proved. ■ 



IV. Optimal Solution for the DC Case 

In this section, we solve Problem (PI) for the DC case. 
We first present a monotonic property for the optimal power 
allocation in Problem (PI), upon which we then develop an 
efficient algorithm to solve this problem. 



A. Monotonic Power Allocation 

Since the optimal solution of Problem (PI) may not be 
unique (cf. Remark 13.11 ), we are interested in finding one 
optimal solution for this problem that leads to the minimum 
power consumption at the source and relay. For such an 
optimal solution, we have the following monotonic property. 

Proposition 4.1: The optimal solution of Problem (PI), 
satisfying ( TTBT l for the case of < ho < 1 or (TT6b for the case 
of ho = 0, is non-decreasing over i, i.e., Pg(i) < Pg(i + 1), 
P R (i+l)<P R (i + 2),i=l,2,--. ,N-l. 

Proof: See Appendix lAl ■ 
With the above monotonic properties, we next address the 
optimal solution of Problem (PI) for the cases with and 
without the direct source-destination link, respectively. 

B. The Case With Direct Link 

For the case of < ho < 1, we consider the Lagrangian of 
Problem (PI), which is given by d28i i. where pk, Afe, ji, and 
are the non-negative Lagrangian multipliers. By taking 
the derivative over Ps(i) and P R (i+l), d29b-d30b are obtained. 
Then, by letting ^^-y = and 9Pj f ( f +1) = 0, the optimal 
solutions for Problem (PI) are obtained as follows: 

(l-fto)Pj(i) 
l+fc P|(i) ' 



1) Case I: if P*(i + 1) > 



p s^ = (i( 



P* (i J- 1 1 — (!-M P s(') 
r R\ l + L )- l+h PS(i) 



- 1 



(30) 



sW 

(l-ho)PJ(Q 



2) Case II: if P*(i + 1) < 1+hoP . {i) 
( 




4(^+1) £™ =I w= 
1 



v+ ■ (3D 

,4(w+i)i:X~^ ~ 1 ) 

Remark 4.1 : From the above expressions, it is observed 
that the source and relay power profiles need to be jointly 
optimized, since the achievable rate for the i-th source message 



l—<£c{h<>Ps(i))+c(p R (i + i)) 

s. t. C(h Q P s (i)) +C(P R (i + l))<C(P 8 (i)),i = !,••• ,N,md ©. 



(PI*) max 

{Psm,{p R (L+in 2(N + i) ^ 



(26) 
(27) 



1 N 

C(Ps(i),PR(i + l),fik, A fc , 7l , Vi+i ) = 2( N + i) £ min {C (Ps ( * } } ' C ( hoPs W ) + C ( Pr (j + X) ) } 



JV 



N 



N 



fe=l \i=l 



fe=l 



dC 



dP s (i) 

dC 

dP R (i + 1) 



4(JV+1) 
1 


X 
X 


1+Ps(») ' 
h 


4(JV+1) 


l + h„P s (i) ' 


o, 






i 


X 


1 



(l-/t )Ps(i) 



B^2^k +7t. 

fc=i 
N 



4(N+1) ~ 1+P fl (i+1)' 



, otherwise 
otherwise 



(29) 
(30) 
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is limited by the available source power in Case I, but by the 
available relay power in Case II. 

From the KKT optimality conditions of Problem (PI), it 
follows that Xk and /ik are strictly positive only when their 
corresponding relay and source energy constraints are satisfied 
with equality. Thus, it follows that the optimal source power 
can change the value from one block to another only when the 
harvested source energy is exhausted at the current block or 
there is a transition between the two values given by ( f30l > 
and ( Bil l. The latter case is due to the fact that ho < 1, 
and thus the source power values given by (130b and (l3"TT l are 
different even when the source energy constraint is not active. 
By further considering the result of Proposition 14.11 we know 
that changing source power values from ( f30b to fl3TT l is not 
possible, and only transitions from (|3TT l to d30T > can occur. 
Thus, we have the following proposition for the optimal source 
power allocation. 

Proposition 4.2: Consider the optimal source power 
P^(i)'s for Problem (PI), which satisfy Proposition 
14.11 For any two successive source energy exhausting 

;., the source 



respectively, which are given as 



ir,o = ar g mm 



mi<j<N { J 



(32) 



with E R (i+l) denoting the relay energy left before the (i+1) 
th relay transmission block, i.e., E R (2) = and E R (i + 1) = 



El=i E R( k + 1) - BP R ( k + 1). * = 2, ■ • • , N, and 



p; 



£.0 



P, 



i+1,0 



(i a , -i + l)B' 



R 



(irfi 



l)B 



(33) 



(34) 



Next, by comparing P s ' a and P R + ' , we determine which 



scenario shown in Proposition 14.21 should happen: 

1) If P* +1 '° > {1 ~ ho) ^ " , it is claimed that Scenario I 
happens, and the optimal source and relay power values are 
given as 



blocks, hi and fc 7 with k% < 
energy constraints 2~2i=i P s(v 



k j ■ 

< 



aR d ^ sE^i^sW 3X6 active, while 

the other constraints Z~2i=i P s(^) < sEL^sf')' 



P*s(n) = P's 
P R (n+l 



i.O 



(l-h )P«' u 



(35) 



l+h P's l 

Then, we set i = i s o + 1, and continue the forward search. 



J 



1, 



1, are all inactive, the optimal source 2) if p i+1 '° < 



power values from the (fcj + l)-th to the kj -th blocks can 
only be one of the following three cases: 



(l-h )P' s ' 
l+hoP^ 



,kj, are identical and 
, kj, are identical and 



1) Scenario I: Pg(i), i — k, + 1. 
given by d3Qb ; 

2) Scenario II: Pg(i), i = h + 1 
given by ( |3T| i: 

3) Scenario III: There exists ko with ki < ko < kj such that 
P|(i) = P , i = h + 1, ko, and P*(i) = P + J- - 
1, i = ko + 1, ...fcj, where Po is jointly determined by 
(l30l > and OTI ). Define the fco-th block as the source power 
transition block for this scenario. 



Scenario II or III may happen. 
To determine whether Scenario III happens or not, we need 
to compute the index ko, 1 < ko < N, of the source power 
transition block defined in Scenario III of Proposition 14.21 

Let the indices i s ,k and i r , p , k,p > 1, correspond to 
the source and relay energy exhausting blocks after i St k~i 
and v )P _i, respectively, and Pg k and P R +1 ' P be the source 
and relay power values between the two successive energy 
exhausting blocks, i.e., 



Based on Proposition l4.2l we know that if we could identify 
all blocks at which the source energy gets exhausted and 
furthermore all the scenarios corresponding to Proposition |4j2] 
the optimal source and relay power profiles for Problem (PI) 
can be obtained accordingly from ( f30b and d3TT i. Thereby, we 
propose Algorithm U summarized in Table U to solve Problem 
(PI), whose optimality proof is given in AppendixIBI The main 
procedure of this proposed algorithm is described as follows. 

Starting from the first block, Algorithm Q] implements a 
forward searching for the optimal power allocation until the 
A^-th block is reached. Suppose that the (i — l)-th block is an 
energy exhausting block for the source, and the optimal power 
allocation for the source and relay have been obtained up to 
this block, denoted by Pg(n) and P R (n + 1), respectively, 
n = 1, •■ • ,i — 1. Then, starting from the i-th source and 
(i + l)-th relay transmission blocks, we first compute i s ,o 
and i r< o (defined in d32l)). corresponding to the next possible 
source and relay energy exhausting blocks, respectively, and 
the source and relay power values Pg° (from the i-th to i s .o-th 
blocks) and P R +1 '° (from the (i + l)-th to (i r . + l)-th blocks), 



i s k = arg mm 

U,k-i<j<N 



i r p = arg mm 

ir, P -l<j<N 



Ci-i.,*-i)B J ' 



(j - i r , p -i) B 



p; 



i . k 



E^ s , fc _ 1+1 ^(fc) 



^ s ,k 



-i)B ' 



?k+1 . p = E^_ 1+1 ^ + D 

[ij-^p ^r.p—l) B 

Then, we find ko > i such that 

5 /• , !^ (1 - ho)Ps(j) 
Pr{3 + 1) < . ~ — , 3=1- 



P R (ko 



2) > 



l + hoPs(j) 
(1 - h )P s {ko + 1) 

l + ^ Ps(fc + l) ' 



(36) 
, (37) 
(38) 
(39) 

(40) 
(41) 



where P s (j) = P s ' k for i s ^ 1 < j < i Sjk , and P R (j + 1) = 



for i r , p -i < j < i r .p, k,p > 0, assuming i s ,-i 



pi+l,p 
ir, — 1 = t — 1. 

Suppose that ko is found in the above. Assume that the 
source power values from block i to fco are P|, and from block 
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fco + 1 to the next source energy exhausting block denoted by 
■here j s is given by 



j s are P s + ^ — 1, where j s is given by 



j s = arg mm 

feo+i<i<JV 



and P l s is given as 



(j -i + l)B 



PL- 



(42) 



(43) 



s {j s -i + l)B 

Then, we further check the following conditions inspired by 



Proposition 14.11 and Scenario III in Proposition 14.2 



Ps ° >H> max ■ 



Ph + T- - 1< 



Pi?(fc + 1) 



1 - h - h P R (k + 1) 
Pfl(fc + 2) 



fto " 1 - h - h P R (k + 2) 

If (l44l and d43T > are true, it is confirmed that Scenario III 
happens, and the optimal source and relay power values are 
given as 



p: 



s- 



,ko 



^s + T^- 1 , ™ = ko + !,••• J. 



(46) 



^r.p—l ^ Tl ^ lr,p> 

PUn + 1) = I n<k ,p>0 . (47) 

l+h„P'(n) ' ' b ^0 T ^ , Js 

Then, we set i = j s + 1, and continue the forward search. 

If k satisfying (|40b and d4TT > cannot be found, or with such 
k Q found in (gUJ and gT} but the conditions in ([44j and (|43T > 
are not satisfied, we claim that Scenario III cannot happen and 
Scenario II must be true. In this case, the optimal source and 
relay power profiles are given as 

5»+i,o 



P*s(n) 



PI' , P R (n + l) 



P 



n = i, ■ ■ ■ ,min(fco,i s ,o), (48) 

where fco = oo if no such ko satisfies ( f4Qb and (PiTl i. Then, we 
set i = min(fco, i Sl o) + 1, and continue the forward search. 

Remark 4.2: According to the proof given in Appendix 151 
it follows that for the case of < ho < 1, the optimal source 
power solution of Problem (PI) obtained using Algorithm Q] 
is unique, while this is not necessarily true for the obtained 
optimal relay power solution according to Proposition 13.11 
However, the obtained relay power solution achieves the min- 
imum energy consumption at the relay, since <fl~5T > is satisfied 
for each i, i = 1, ■ ■ ■ , N. 

C. The Case Without Direct Link 

Similar to the case with direct link, we obtain the following 
optimal power solutions for Problem (PI) in the case of ho = 
0: 

1) IfP*(i + l)>PJ(i), 



p* R (i+l) = p* s (l) 



- 1 



(49) 



TABLE I 

Algorithm[2 Compute the optimal solution of Problem (P 1 ) 

FOR THE CASE OF < ho < 1. 

• Set i = l; while i < N, repeat 

• Compute i Si o and t rj o by )32t . and Pg° and P^ 1 ' by (33) and (34). 

1) If P'j+ lfi > (1 ~ ho) ~S ° , compute P*(i) and P* (i + 1) by (35), 
and set i = i s> rj + 1; 

2) If P^ 1 ' < ^ — fc°^ S ° ' com P ute *s,fe> V,p. Ps' k ar, d P^ 1,p , 
k,p>l by (36)-(^), S and check 

a) If there exists fco satisfying )40t and )41t , compute j s and Pg 
by (42) and (43), respectively, and check 

i) If P^ satisfies (44) and (45), compute Pg(i) and P^(i+l) 
by (46) and (47); set i = j 3 + 1. 

ii) else compute P£{i) and P^(i + 1) by (48) ; set j = 
min{fco, i s ,o} + 1. 

b) else compute Pg(i) and P^(« + 1) by (48); set j = i Sj o + I. 

• Algorithm ends. 



2) If P R (i+l)<PI(z), 

r p|(i)=^+i) 

l ^tt(i + l) = 



1 



(50) 



It is worth noting that to obtain the optimal power profile with 
the minimum energy consumption, we set the source and relay 
power levels the same, while in general this is not necessary 
since the optimal source and relay power profiles may not be 
unique. 

From the above solutions, it is observed that: 1) The source 
power Pg(i) at the i-th block and relay power P R (i + l) at the 
(i + l)-th block should be identical, i = 1, ■ • ■ , N; 2) Since 
Xk and fik are strictly positive only when their corresponding 
energy constraints are satisfied with equality, it follows that 
the source/relay power changes value only when the harvested 
energy at either the source or the relay is exhausted. 

Based on the above observations, Algorithm U for the case 
with direct link can be simplified to obtain the optimal source 
and relay power allocation for Problem (PI) in the case 
without direct link. We denote this algorithm as Algorithm HI1 
which is summarized in Table [TTJ Since the optimality proof 
of Algorithm II is similar to that of Algorithm I, it is omitted 
here for brevity. 

Remark 4.3: The optimal solution for Problem (PI) ob- 
tained using Algorithm [TT] in the case of ho — is illustrated 
in Fig. |2] where we see that the optimal source/relay power 
profile corresponds to the shortest path from the origin to 
the highest possible sum-energy point at the end of the A^- 
block transmission under two stair-like source and relay energy 
constraints. It is worth noting that for the case of ho — 0, 
our relay channel model degrades to a cascade of two AWGN 
point-to-point channels that were studied in Q with individual 
EH constraints. Thus, at each block, either the source-relay or 
the relay-destination link can be the performance bottleneck. 

V. Optimal Solution for the NDC Case 

In this section, we solve Problem (P2) for the NDC case. 
We first prove that a separation principle for the source and 
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TABLE II 

AlgorithmHH Compute the optimal solution of Problem (PI) 

FOR THE CASE OF ho = 0. 



• Set i = l; when i < N, repeat 

• Compute Pi,, P^ 1 , is, and i r as follows. 



i s = arg mm 

i<j<N 



i T = arg mm 

i<j<N 



<J ~ i + 1) B J ' 



(j-i+l)B 



s {i s - i + 1) B 

fl (i r - i + 1) B 

where B s (l) = E R {2) = 0, = T,t~Jl E s(k) " SP^(fc) 

and B K (t + 1) = El=i + 1) - BP* (fc + 1), i = 2, ■ ■ • , N, 

respectively. Then, the optimal source and relay power profiles are given 
as 

1) If P» > P^ +1 ,setP*(fc) = P R {k+l) = P^ +1 for k = i, ■ ■ ■ ,t r , 
and_i = + 1; 

2) If P| < PJ+ 1 , set P|(jfc) = P^(fc + 1) = P| for it = i, ■ ■ ■ , i s , 
and i = i s + 1, 

Algorithm ends. 



Source energy profile 
Relay energy profile 
Optimal power profile 

PI 




Fig. 2. An example for the optimal source and relay power allocation for the 
DC case with ho = 0. Note that the slope of each dashed line corresponds 
to Pg or P^~ obtained using Algorithm [n] 



relay power allocation problem holds, upon which Problem 
(P2) can be solved by a two-stage strategy: First obtain the 
optimal source power allocation by ignoring the relay, and then 
optimize the relay power allocation with the obtained source 
power solution. Since this separation principle applies to both 
cases with and without the direct link, we address these two 
cases with a unified analysis as follows. 

A. Optimal Source Power Allocation 

First, we consider the following source power allocation 
problem by ignoring the relay: 



A' 



(P3) max y^C(hPs(i)) 



(51) 



s.t. J2 P sW ^ ^E^W, k= 1,2, • •• ,N, (52) 



where h is a constant with < h < 1. Problem (P3) has been 
solved in 0, for which the algorithm (denoted by Algorithm 
HIIl l is summarized in Table [III] for the sake of completeness. 
Note that the optimal source power profile PS(i)'s of Problem 
(P3) are non-decreasing over i |3). 

table m 

AlgorithmHTII Compute the optimal solution for Problem (P3). 



• Initialize i = 1; while i < N, repeat 

• Compute 

. [ J2U t E s {k) \ J2lu E s(k) 

i s = arg mm < , ^ — \ PI = - — E=2 — . 

i<i<N\(j -i + l)B j ' s {i s -i + l)B 

The optimal source power profile is given as Pg{n) = Pg, n 
i, ■ ■ ■ , i s ; set i = i s + 1. 

• Algorithm ends. 



It is easy to see that the optimal source power profile of 
Problem (P3) maximizes the average throughput of both the 
source-relay and source-destination links. Moreover, since for 
the NDC case, the relay can store the binning indices of the 
decoded source messages with arbitrary delay before forward- 
ing them to the destination with best effort transmissions, the 
relay power profile intuitively should have no effect on the 
optimal source power profile. This conjecture is affirmed by 
the following proposition. 

Proposition 5.1: For the NDC case, the optimal source 
power solution for Problem (P3) is also globally optimal for 
Problem (P2). 

Proof: See Appendix lei ■ 
This proposition implies that the separation principle for 
the source and relay power allocation problems is optimal for 
Problem (P2). Thus, even though Problem (P2) is non-convex, 
we can still find its globally optimal solution efficiently. 

Remark 5.1: It is shown in Appendix ICl that for the case 
of < ho < 1, the optimal source power profile Pg(i)'s of 
Problem (P2) given by Algorithm [HI] is unique. However, for 
the case of ho = 0, this result is not true in general, since the 
source energy may not necessarily be exhausted at the end of 
each iV-block transmission for Problem (P2) (cf. Proposition 

B. Optimal Relay Power Allocation 

With the optimal source power profile P|(i)'s obtained 
using Algorithm [Hi] the optimal relay power profile can be 
determined by the following problem: 

A 

(P4) max VC(Pr(;+1)) (53) 

P H (i+l)>0,Vi ^ 

1=1 

k k 

s.t. ^C(Ph(* + 1))<E C ( P s«) 



! = 1 



i=l 



J2C(hoPs(i)), k = l,---,N, and,©. (54) 



This problem is non-convex due to the constraints in d54"l ). 
However, by letting r(i + 1) = C (Pr(i + 1)), Problem (P4) 
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can be rewritten as 



Rate 



iY 



(P5) max VVfi + l] 

r(i+l)>0,Vi ^ 



(55) 



s. t. J2r(i + l)<J2C(P£(t))-C(h P* s (i)), (56) 

i=l i=l 

E(2 2r(i+1) -l)< ^^(i + l), 



k = l,---,N. (57) 



It can be shown that Problem (P5) is convex over r(i + l)'s. By 
the KKT optimality conditions, we obtain the optimal solution 
for Problem (P5) as 



1 , 1 - Efc=j X k 



r*{i + l)= -log 



2 21n2.££V7 fe 



(58) 



where Aj. and 7^ are the non-negative Lagrangian multipli- 
ers corresponding to the fc-th constraint in (l56l l and (fSTb . 
respectively. Problem (P5) can be solved by a forward search 
algorithm, denoted by Algorithm [IV] in Table |IV] for which 
the optimality proof is similar to that of Algorithm HU and is 
thus omitted. 

TABLE IV 

AlgorithmUVI Compute the optimal solution for Problem (P5). 



Initialize i = 1; while i < N, repeat 
Compute 



1 1 = arg mm 

i<j<N 



12 = arg mm 



Ci + J2i^c(p* s (k))-c(h p*(k)) 

(j _ i + i) B 



E i+1 + -£l =i E s (k + l) 



ri 



<j<N (j - i + 1)B 



r 2 =C 



(h-i + l)B 
(i 2 -i + l)B 



where Ci = E 2 = 0,0,= J2' k ~Ji C (P£(k)) - C (h PS(k)) 

r*(k), and E i+1 = Efc=\ (s S (* + 1) - 2 2 ''*( fc + 1 ) - l), i 

2, • • ■ , AT. Let jo = argmin J= i,2 {rj}. Set r*(j + 1) = r jo , j 

»."• »»j - and * = »jo + L 
Algorithm ends. 



Remark 5.2: It is worth noting that from (l58l l. we observe 
that the optimal relay transmission rate r*(i + 1) is non- 
decreasing over i, and strictly increases when any one of 
the constraints (l56l ) and (l57l i is satisfied with equality. Thus, 
the optimal relay power profile P^(i + l)'s of Problem (P2) 
with P£(i + 1) = 2 2r *< l+1 ) - 1, i = l,-- - ,N, are also 
non-decreasing over i. Furthermore, it can be shown that the 
obtained P^(i+l)'s achieve the minimum energy consumption 
at the relay. 

Based on the analysis in the above two subsections, we 
obtain the following proposition for the monotonic properties 



Y/Z/Z/A _ C(/- ( ( )) C(M'(„) 



J 1 N 

Fig. 3. Illustration of the rate allocation for + l)'s in the NDC case. 



for Ps(i)'s and P^(i + l)'s for the NDC case, similar to the 
DC case (cf. Proposition 14. It . 

Proposition 5.2: For Problem (P2), there exist optimal 
source power solution Pg (i)'s and relay power solution P^(i+ 
l)'s that are non-decreasing over i, i = 1, • • • , N. 

Remark 5.3: It is worth noting that problems of similar 
structures as (P4) have been solved independently in prior 
works [3], where the results are based on generalizing the 
solution for the two-epoch case. Nevertheless, in this paper, 
we use some optimization tools to transfer Problem (P4) into 
Problem (P5), which is convex and much easier to be solved. 



C. Optimal Rate Scheduling 

It has been shown in Section IV that for the DC case, the 
relay binning rate Rb{i + 1) for each source message can 
be directly computed by ( fTTI i given the optimal source and 
relay power profiles. However, for the NDC case, it requires 
additional effort to obtain Rs(i+1), since we need to consider 
two sets of different constraints in (f20t and (1221 (as will be 
shown next). Suppose that + 1) is obtained with the 

optimal source and relay power profiles Pg(i)'s and P^ii + 
l)'s, the source transmission rates R(i)'s can be determined 
from (fJTJ. Then, with the obtained R(i)'s, we can update the 
obtained source power profile to achieve the same maximum 
throughput but with the minimum energy consumption. 

To compute + l)' s > the following observations are 

first drawn. If C(P^(i + 1)) > C(P£(i)) - C(h P*(i)), Vi £ 
{1, • - • , N}, the relay should transmit not only the binning 
index of the i-th source message at the (i + l)-th block, but 
also those of source messages 1 < j < i. Moreover, due 
to the constraint in (f25]l, it follows that if C(P^(i + 1)) > 
C(Pg(i)) — C(hoPjs(i)), Vi, there must exist j with 1 < j < i, 
such thatC(P£0' + l)) < C(Pg(j))-C{hoP£(j)). The above 
observations imply that to obtain Rb{i + l)'s, we need to 
find all i's with C{P^{i + 1)) > C(P|(i)) - C(h P£(i)), and 
then use their surplus rates to transmit the binning indices 
of source messages j < i. In Fig. [3] we show an example 
for the relationship between C(P^(i + 1)) and C(PJ(i)) - 
C(h Pg(i)), where Case A stands for the case of C(P^(i + 
1)) < C(P^.(i)) - C(h Pg(i)) and Case B for the case with 
a reversed inequality. Then, a geometric interpretation for the 
algorithm of computing Rb{i + l)'s is given as follows: Use 
the surplus rate area in Case B to fill the deficient rate area in 
Case A. Obviously, if the area to be filled in Case A is larger 
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than the filling area in Case B, the corresponding values for 
Rb{i + l)*s are not unique. 

Thus, we now develop a backward search algorithm (Al- 
gorithm |Vj that is summarized in Table [V] to obtain one of 
the feasible solutions for Rs{i + l)'s. The main procedure 
of this algorithm is described as follows. First, Rb{i + l)'s 
are initialized as the minimum values between C(P^(z +1)) 
and C(Pg(i)) — C(hoPg(i)) for all i's, and a parameter t 
(which, e.g., corresponds to the filling area in Case B of Fig. 
3) is set to be 0. The algorithm then searches the values for 
Rs{i + l)'s in a backward way from i = N to 1. For any 
i-th block, the algorithm computes temp = C(Pj£(i + 1)) — 
(C(PJ(i)) - C(Ms(i))). Then, if temp > 0, Case B occurs, 
and temp is added to t; if temp < 0, Case A occurs, and 
Rb(i + 1) is increased by min(— temp, t), and this amount is 
then subtracted from t. 

TABLE V 

AlgorithmIvI Compute R b (i + l)'s for the NDC case. 



. Given Pt(i)'s by Algorithm |m] and P^(i + l)'s 
by Algorithm [TV] initialize fi B (i + 1) = 
min {C (P£(i + 1)) ,C(P|(i)) -C(hoP$(i))}, i = 1,---,N, 
and t = 0; 

• From i = N to 1, compute temp = C (P^(i + 1)) - C(Pg(i)) + 

C{hoP§{i)) 

1) If temp > 0, set t = t + temp; 

2) If temp < 0, set Rb (i) = Rb (*) + min(— temp, t) and t = 
(t + temp) + . 

• Algorithm ends. 



Remark 5.4: It is recalled in Remark 15.11 that in the case 
of ho = 0, the optimal source power profile Pg(i)'s given by 
Algorithm Hn] may not achieve the minimum energy consump- 
tion. In this case, from ( f2TT >. it follows that R(i) = Rb{i + 1), 
i = 1, • • • ,N, with Rb{i + l)'s obtained by using Algorithm 
|V1 In order to achieve the minimum energy consumption at 
the source, the optimal source power solution of Problem 
(P2) can be updated as P|(i) = 2 2Rb{i+1) - 1. It is worth 
noting that the above obtained source power profile is still 
non-decreasing over i (thus in accordance with Proposition 
I5.21 i. since Rg(i + l)'s obtained using Algorithm [Vl are non- 
decreasing over i. 



D. Throughput Comparison: DC vs. NDC 



As shown by Proposition 13.21 the throughput of the NDC 
case is no smaller than that of the DC case. To further compare 
these two cases, the following proposition shows when the 
NDC case is strictly better than the DC case in terms of the 
average throughput. 

Proposition 5.3: The average throughput of the NDC case 
is strictly larger than that of the DC case if and only if there 
exists i g {1, • • • , N} such that C(P^(i + 1)) > C (P|(i)) - 
C (hoPg(i)), where P|(i)'s and C(P^(i + l))'s are the optimal 
solutions to Problems (P3) and (P5) obtained by Algorithms 
IITI1 and HVl respectively. 

Proof: See Appendix ID1 ■ 
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Fig. 4. Throughput comparison of various power allocation schemes for the 
relay channel with energy harvesting constraints. 



VI. Numerical Results 

In this section, we present some numerical results to validate 
our theoretical results. For the purpose of exposition, we 
assume a periodic energy profile model for some predictable 
EH sources. Specifically, the source and relay energy profiles 
are given as 



( ^T 2?r + f 



Ac 



A, 



1< i < N, 



respectively, where As,Ar > are the amplitudes of the 
sinusoidal energy profiles at the source and relay, respectively, 
and 6 is the phase shift between these two energy profiles. 
Here, we choose B = 100, N = 40, 6 = §tt, and A s = 
A R = 200. 

We compare our proposed algorithms with a greedy power 
allocation strategy. Here, we adopt a non-trivial greedy algo- 
rithm by assuming that both the source and relay know the 
harvested energy amounts up to the current block prior to 
transmissions. The transmission rate of the i-th source message 
in the greedy algorithm is given as 




(59) 



where E s (i) = £Li E s (k) - PEl=i Ps{k), and E R {i + 
1) = ELi E R( k + 1) - B El~Ji Mb + 1), with P s (i) = 



2 R G (i)-C[ h 



Es(') 



2 2fl G (i) _ i and p R (i + i) = 2 v ~ r ° 

i — 1, • • • ,N. Note that in the above greedy algorithm, both 
the source and relay consume as much available power as pos- 
sible at two successive blocks to maximize the instantaneous 
throughput, thus achieving the minimum delay as for the DC 
case. 

In Fig. |U we show the average throughputs versus the 
direct link channel gain ho for the proposed power allocation 
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algorithms and the greedy algorithm. It is observed that as 
the direct link becomes stronger, i.e., ho increases, there is 
a throughput limit of 0.387 bps/Hz. For the NDC case, this 
throughput limit is achieved even for very small ho around 
0.05. In contrast, for the DC case, the throughput increases 
almost linearly and achieves the throughput limit when ho 
exceeds 0.75. Thus, the throughput gain of NDC over DC 
cases by relaxing the decoding delay is present only in the 
regime of small direct link gains, thanks to the exploitation of 
energy diversity in cooperative communication. Also note that 
with small ho, the condition in Proposition 15.31 is more likely 
to be satisfied. For the greedy algorithm, it is observed that 
the throughput loss can be large, especially when ho is small, 
as compared to the proposed algorithm for the DC case. 

VII. Concluding Remarks 

In this paper, we studied the throughput maximization 
problem for the orthogonal relay channel with EH source and 
relay nodes, assuming a deterministic EH model. For both the 
cases with and without delay constraints at the destination, 
we examined the structures of the optimal source and relay 
power profiles over time, and developed algorithms to effi- 
ciently compute these optimal power profiles. In addition, a 
new interesting energy diversity phenomenon was explored in 
EH-powered wireless cooperative communication with delay 
tolerance. We now conclude the paper by highlighting some 
important aspects unaddressed yet and thus worth being in- 
vestigated in the future work as follows: 

1) This paper considers the deterministic EH model, while in 
practice, many EH sources should be modeled as random 
processes, e.g., the block-Markov model |6j. The study 
of the throughput maximization problem for the Gaussian 
relay channel under random EH models motivated by 
the results obtained here will thus be practically more 
appealing. 

2) To simplify the analysis, it is assumed in this paper that 
the energy storage capacity is infinite, which may not 
be true in practice J4)- Thus, considering the case with 
finite energy storage will be an interesting extension of 
this paper. 

3) In this paper, the DF relaying technique is adopted to 
design the optimal transmission for the orthogonal relay 
channel. Extensions to more general relay channel models 
and/or other relaying techniques, such as amplify and 
forward (AF) as well as CF, will be also interesting. 

Appendix A 
Proof of Propositions. II 



Denote the optimal source and relay power profiles of 
Problem (PI) as Pg(i)'s and P R (i + l)'s, respectively. For any 
consecutive power pair, consider the following three cases: 

1) P|(i) > P£(i + 1) and P*(i + 1) > P R {i + 2): Define 
a new power allocation profile as Ps{i) = Ps('i + 1) = 
p »«+p+ 1 > and P R (i + l) = P R (i+2) = afi±a+Sg±31 , 
It is easy to check that the new power allocation profile still 
satisfies the energy constraint in (fT~4T >; moreover, since the 
objective function in (fT3l ) is concave, it follows that the new 



allocation leads to a larger sum rate over these two blocks. 
Thus, this case cannot happen. 

2) PJ(i) > P*(i + 1) and P R (i + 1) < P R (i + 2): First, 
we prove that the i-th constraint in ([TBT l is not achieved with 
equality by contradiction as follows. Assuming that its equality 
is achieved, it is then observed that P R (i + 1) < P R (i + 2) is 
contradicted with the following results: 

1 + Pg(») 



log(l + P£(i + l))=log 



l + h Q P*(i) 



>l0g|7 i±sa±iL, (60) 



i + h p*(i + i) 

>log(l + Pj^ + 2)), (61) 

where doOb is due to the fact that log ( \+h o3 ^ j is strictly 
increasing over x > for any fixed ho with < ho < 1, 
and doTT i is by Proposition 13. II 



Thus, there must exist 5, < 5 < 



Ps(0-P|(»+1) 



, such that 



C(Pt(i)-6)>C{ho(P£(i)-5)) + C(P R (i + l)). (62) 



Define a new power allocation profile as Ps(i) = Pg(i) — S, 
P s (i + 1) = P|(i + 1) + 6, P R (i + 1) = P R (i + 1), and 
P R (i + 2) = P R (i + 2). It is easy to check that the new 
allocation satisfies the energy constraints in ( TT4l i. and yields a 
larger sum rate over these two blocks, i.e., 

R(i) + R(i + 1) 

> C (ho? s (ij) + C (P R (i + 1)) + C (h P s (i + 1)) 

+ C(P R (i + 2)) (63) 

> C (h P£(i)) + C (P R (i + 1)) + C (h P£(i + 1)) 

+ C(Pj^ + 2)) (64) 
= R*(i) + R*(i + 1), 

where d63l is due to (l62l and the following fact: Since C(x + 

6) -C(x) > C (h (x + 6)) - C (h x) for < h < 1, it 
follows that 



C {P 8 {i + 1)J - [C [hoPs(i + l))+C (P* R (i + 2))^ 

>C (P|(i + 1)) - C {hoP* s {i + l))-C (P R (i + 2)) > 0; 

and (164-b is due to the concavity of log(x) and 5 < 
£s_C0 J's('+ 1 ) _ Therefore, this case cannot happen. 

3) P£(i) < P|(i + 1) and P*(i + 1) > P R (i + 2): Since 
l°g f i+feox ) i s strictly increasing over x > for fixed ho, 
< ho < 1, it follows that 



lop 



1 + h Q P*(i + I) J ~ + h a P*(i) 

>log(l + P£(i + l)), 



where the last inequality is due to Proposition 13.11 Thus, it is 
obtained that 

C (P* s (i + 1)) > C (h P* s (z + l))+C (P R (i + 1)) 

>C(h P$(i + l)) + c(P R (i + 2j) , (65) 

where P R (i + l) = P R {i + 2) = P ^ (i+1) + P « ( ' +2) . By keeping 
Pg{i) and PZ(i + 1) unchanged and updating the relay power 
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values as P R (i + 1) and P R (i + 2), it is observed that the 
relay energy constraints are still satisfied, and the sum rate is 
improved, i.e., 

R(i) + R(i + 1) (66) 

= C(h P$(i))+c(p R (i + l) y ) 

+ C(h P£(i + l))+c(p R (i + 2)) (67) 

>C(h P^)+C(P R (i + l)) 

+ C{h P*{i + l))+C(P R (i + 2)) (68) 
= R*(i) + R*(i + l), 

where (|67]) is due to (|65]) and the fact that C(Pg(i)) > 
C (h P£(i))+C (P R (i + 1)) > C {h P$(i))+C (p R (i + 1)), 
and d68l is due to the concavity of log(a;). Thus, this case 
cannot happen. 

To summarize, since all the above three cases cannot be true, 
the only remaining case of Pg(i) < P|(i+1) and P R (i+l) < 
P R (i + 2) must be true. Proposition 14. II is thus proved. 

Appendix B 
The Optimality Proof of Algorithm Q] 

First, we prove that the source power profile Pg(i)'s 
obtained using Algorithm Q] are optimal for Problem (PI). 
Given Pg(i)'s, assume that the equalities of the source energy 
constraints are achieved at blocks i\, is, - • • ,i m = N, while 
those are not achieved for the other blocks. Moreover, we 
assume that before block i s , < s < m (define %q = 0), 
the optimal solution Pg(i)'s of Problem (PI) are the same as 
Pg(i)'s, and their difference first appears at the ?'-th block, 
i s < i < i 8 +i. Due to this difference, the index of the next 
source energy exhausting block is denoted as i s +i, which may 
not be equal to i s +i- Then, only three scenarios shown in 
Proposition 14.21 may happen for both P^(z)'s and P^(i)'s for 
the (i s + l)-th to the i s+ i-th blocks and the i s+ i-th block, 
respectively, which are discussed as follows: 

(I) Scenario I happens for P^(i)'s, and consider the fol- 
lowing two cases: 

(a) If Pg(i) > Pg(«) : B Y propositionl4~n it follows that 
P|(j) > P|W > P|(<) = P|(i), 3 = i, ■ ■ ■ , 

and then B^J^W > B PJ (*) = 
Y^ik=l Es(k), which violates the source energy 
constraint. Thus, this case cannot occur. 

(b) If P|(i) < P£(i): For the case that Scenario I or II 
happens for Pg(i)'s, it is easy to check that there 
will be no block i s+ i > i where the source energy is 
exhausted, and thus Pg(i)'s violate Proposition [33] 
For the case that Scenario III happens for P|(z)'s, it 
is observed that both the fco-th source and relay en- 
ergy constraints in (Q]i and (f2|i are not achieved with 
equality, where fco is the index of the source power 
transition block shown in Scenario III of Proposition 
14.21 Then, there must exist < 9 < 1, such that the 
new source and relay power allocation profiles de- 
fined below still satisfy the source and relay energy 
constraints: P s {k ) = 6P£(fc ) + (l-6>)P^(fc + l), 



Psikp + 1) = (1 - O)PI(ko) + 9P*{k + 1), 
Pfl(fco + l) = eP R (k + l) + (l-e)P R (k + 2), and 
P R {k Q + 2) = (1 - 9)P R (k + 1) + 9P R (k + 2). It 
is easy to check that with the new power allocation 
profiles, the sum rate over these two blocks is 
improved, since the rate function dTZb is concave. 
Thus, this case cannot occur. 
(II) Scenario III happens for P|(i)'s, and assume that 
PsU) = Po, J = is + 1, •••,*(>, and P*(j) = 

P + ^ 1, j = fco + , »g+i- Consider the 

following two cases: 

(a) If Pg(£) > P|(«) : When i > i s , where i s>0 is 
given by d32l . it is easy to check that only Scenario 
III can happen for PS(i)'s from the (i s + l)-th 
to the i s+ i-th blocks. However, the source energy 
constraint at the i s+1 -th block is violated, since 
Pg(z)'s are non-decreasing by Proposition 14.11 As 
such, we only consider the case of i < i s , which 
consists of two subcases: 

(1) P|(i) < P£(i) < P l s :°, where P|° is given 
by ( f33l > and d34b : For Pg(i)'s, similar to case (lb), 
it follows that only Scenario III can happen, and 
denote fco as the index of the source power transition 
block defined in Scenario III of Proposition 14.21 If 
fco < fco, by a similar argument as case (la), the 
source energy constraint at the i s+ i-th block will 
be violated; if fco > fco, by a similar argument of 
case (lb), it follows that Pg(i)'s cannot be optimal. 

(2) Pg(i) = Pg° : First, it is claimed that from 
the (i s o + l)-th to the i s+ i-th blocks, there is no 
such index fco corresponding to the source power 
transition block for Scenario III in Proposition 14.21 
This is proved by contradiction, and consider the 
cases of fco > fco and fc < fco following the same 
argument as case (Hal), respectively. 

As such, for Pg(i)'s, it is obtained that Scenario II 
happens from the (£ s + l)-th to the fco-th blocks, and 
Scenario I happens from the (fco + l)-th to the i s +i- 
th blocks. Next, we prove that Pg(i)'s are strictly 
sub-optimal over these blocks. Define a new source 

power profile as P s (j) = — (£-i,) B — ' 3 = 

». + l t - ,k andP g (j) = £ ^ff ) , 3 = 
fco + 1, • • • , i s +i- Since both log(x) and log(h x) 
are concave and some of the source energy con- 
straints at the (fco + l)-th to the i s+1 -th blocks may 
be violated by Ps(j)'s, the sum rate Y^j=i s +i B(i) 
given by Pg(i)'s is upper-bounded by 

E R M< E C(h P s (j))+C(P R (j + l)) 

j=i 3 +l j=i a +l 

is+i 

+ Yl C (Ps(j)), (69) 

j=k + l 

where P R (J + 1) = P R + ,p , with i r , P -i < j < i r ,p, 
P > 0, i r .-i = i s , j < fco> and i rtP 's and P R + ' p 's 
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are given by (O, (03), (O, (O, t®, and <[39]i. 
On the other hand, the sum rate for the (i s + l)-th 
to i s +i-th source messages with P^(z)'s is 

R*= £ C(Aoi^(j))+C(PA(j + l)) 

j=i 3 +l j=is+l 

+ E c ( p sa)), (70) 

3=k + l 

where P£(j + 1) = P^ +1 - p = P R (j + 1) is due 
to (1471 1. To prove that ( l70l i is larger than the right 
hand side of d69l ), it is equivalent to show that ( fTTT i 
is true, where 5 = Ps(i s + 1) - P§(is + 1) = 
^sC^ + l)-^, ii = fco— and Z 2 = i s +i-fc . By 
removing the log operations, (TTTt can be rewritten 
as 



(1 + fto (Po + 5)) h P + i - k<5 



> 1. 



(72) 



Let = (1 + /i (Po+x))' 1 

Note that for < x < S, P + 
thus it follows that 

g'(x) = -h(l + h (P + x)) 



Pi 



h-1 



h 



x > 0, and 



1 



h h 



+ 1 ) h a x < 0, 

'2 



which means that g(x) is decreasing over < x < 
S, and it follows that ( F72b is true. Then, P^(i)'s 
cannot be optimal for this case, 
(b) If P|(i) < P*(i): First, note that if i > i s + 1, 
Pg(z)'s violate Proposition 14.11 Thus, we assume 
that i = i s + 1. Then, for Pg(z)'s, if Scenario I or 
II happens, it is easy to check that the source energy 
cannot be completely consumed at the end of each 
N -block transmission, which violates Proposition 
13.31 Consider now the case that Scenario III hap- 
pens. Denote the index of the source power transi- 
tion block defined in Scenario III of Proposition [472] 
as fco and there are then two subcases: (i) If ko > fco, 
it is easy to check that this case violates Proposition 
13.31 since there will be no source energy exhausting 
blocks after the fco-th block; (ii) if fco < fco, it is easy 
to check that the relay power constraint at the fco-th 
block is achieved with equality with a similar argu- 
ment as case (lb) under Scenario III. Moreover, it 
can be shown that Pg(J) > PgO'), fco < j < fco (if 
not, there will be no source power energy exhausting 
block existed after the fco-th block). Then, it follows 



that p*(j + 1) > ^H^f > ( ;tM^^ > 

Pp\{j + 1) for fco < j < fco- As such, it is observed 
that the the relay power constraint at the fc -th block 
is violated. Therefore, this case cannot occur. 
(Ill) Scenario II happens for P^(i)'s, and consider the fol- 
lowing two cases: 

(a) If P|(i) > P|(i): Similar to case (la), it follows 
that P|(i)'s cannot be optimal. 

(b) If P£(i) < P|(i): For P|(i)'s, Scenario I and II 
cannot occur due to the same argument as case (lb); 
for the case that Scenario III happens for Pg(i)'s, 
it cannot happen, since by using Algorithm U we 
cannot find such fc and Pg satisfying (1401 . fill . 
(1441 . and (1451) . respectively. Then, we only need to 
show that (l40b . WH , (1441 1. and (l45T l are necessary for 
the existence of Scenario III. By a similar argument 
as case (II), it can be shown that (l40l and (HTt are 
necessary for the existence of fco. For (PPfl i and (1451 . 
it is shown as follows: (i) P l s < Py°: this is due 
to the source energy constraint at the i Sj o-th block; 
(ii) Ps > -fj?(*s): this is due to Proposition |4~T1 (iii) 

which is equivalent to that 



pz > Pi? (fco + l) 

S - 1-ho-haPnXko + l) 



P R (k + l) < 



If this condition is not true, 



(i-Mfs . 

1+hoPi. 

which means that the fco-th source and relay energy 
constraints in (Q]i and (f2]i are not achieved with 
equality, it can be shown that this case is not optimal 
with the same argument of ca^se (lb) under Scenario 
III; (iv) P| + i - 1 < 



Pn(fc +2) 



, which is 

l-ho-hoP R (ko+2) 
~ <l-h )(pi + -l) 

equivalent to that P B (k + 2) > '\ h ° J- : 

1+^0(^+^-1) 

If this condition is not true, it will violate the 
optimality conditions for Scenario III shown in 
Proposition 14.21 Therefore, the necessity of (140b . 
fiTJ, (04), and §3]l is proved. 
In conclusion, Pg(iy& are optimal for Problem (PI). Now, 
with the optimal Pg(i)'s, we know which scenario in Propo- 
sition |4j2] happens for each block. It is thus easy to show that 
the corresponding relay power profile PJj(i + l)'s obtained 
using Algorithm U are also optimal: If PJj(i + l)'s are given 
by POl l. they are optimal obviously; if Pjj(i + l)'s are given 
by d3~Tl i. they are also optimal by proving that the two cases 
P£(i + 1) > P* R {i + 1) and P R {i + 1) < P a (i + 1) both 
cannot occur, for which the proof is similar to that in case (I) 
and thus omitted. Therefore, the optimality of Algorithm U is 
proved. 

Appendix C 
Proof of Proposition ^. II 

Define f(x) = \ log (rrfe) over x > with fixed < 



ho < 1. Since f"(x) 



5(1 + 2;) + \h\ (1 + h x)~ 



\C (h Po) + hC Po 



hC(h (P + 6)) + l 2 C P + 



>0, 



(71) 



14 



and 1H ^° - < j^, it follows that f"(x) < and f(x) is 
concave. Moreover, it is easy to check that f(x) is increasing 
over x > 0. Then, we obtain the following lemma. 

Lemma C.l: For the case of < ho < 1, the optimal source 
power profile PJ (i)'s of Problem (P2) are non-decreasing over 
i; for ho = 0, there exist optimal source power profile P|(i)'s, 
which are non-decreasing over i, i = 1, • • • , N. 

Proof: Denote the optimal solution of Problem (P2) as 
P|(i)'s and P£(i + l)'s. For the case of < ho < 1, consider 
any consecutive source/relay power pair corresponding to the 
i-th and the (£ + l)-th source messages, i = 1, • • ■ , N — 1, and 
the following two cases: 

1) P|(i) > P|(i + 1), P£(i + 1) > P%(i + 2): Define 
r*(i + l) = C(Pfj(i + l)). Since log(x) is increasing, it follows 
that r*(i + 1) > r*(i + 2). Then, the constraint in d25l l is 
equivalent to El=i r *(k + 1) < El=i f(Ps( k ))> which is 
convex. Thus, similar to case 1) of Appendix lAl this case is 
not optimal for Problem (P2). 

2) P$(i) > P$(i + 1), P*(£ + 1) < P*(i + 2): It 
is first proved that the constraint Efc=i ^(Pr(^ + 1)) — 
El=i f(Ps(k)) i s not satisfied with equality in this case by 
contradiction as follows. Suppose that the above rate inequality 
is satisfied with equality. From the (i — l)-th and the (i + l)-th 
constraints in ([25), it follows that C(P£(i + 1)) > /(P|(i)) 
and C(P£(i + 2)) < /(P|(i + 1)), and together with the 
assumption that P^(i + 1) < P^(i + 2), it follows that 

< C(F£(t + 1)) < C(P£(i + 2)) < /(P|(i + 
1)). Since f(x) is an increasing function, it follows that 
PgW < P?(* + 1)' which contradicts the assumption that 
P*(i)>P*(i+l). 

Then, there must exist 5 > 0, with which we define a new 
power allocation as Ps(i) = PJ(i)— <5 and Ps(i+1) = Pg(i + 
l) + S, while keeping P^(i + 1) and P£(i + 2) unchanged. It is 
easy to check that the new power allocation still satisfies the 
constraints in Problem (P3) and increases the sum rate over 
these two blocks. Thus, case 2) is not optimal for Problem 
(P2). 

In conclusion, for the case of < ho < 1, only the case 
P|(«) < Ps(« + 1) and P£(i + 1) < P£(£ + 2) can be optimal. 
For the case of ho = 0, similar argument can be applied to 
show the existence of such a non-decreasing optimal solution, 
which is omitted for brevity. Thus, this lemma is proved. ■ 

Next, we prove Proposition 15.11 First, consider the case of 
< ho < 1. For the optimal source power Pg(i)'s of Problem 
(P2) and Pg{i)'s obtained using Algorithm Hill we assume 
that P* s (j) = P|(i), j = 1, • • • ,i - 1, and P£(i) + P|(i). 
Consider the following two cases: 

1) Pg(i) > Pg(i): Similar to the proof of case (la) in 
Appendix iBl it follows that this case violates the source energy 
constraints, since Pg(i)'s are non-decreasing, and thus cannot 
be true. 

2) Pg(i) < Pg(i)' First, we claim the following two 
results, which will be proved later: (i) Pg(i)'s reduce the 
optimal value of Problem (P2) without considering the relay 
power allocation; (ii) By further considering the relay power 
allocation and the constraint d25l ), Pg(i)'s shrink the feasible 
set of Pr(£ + l)'s. If (i) and (ii) are both true, it follows that 
Pg(i) < Pg(«) cannot be true. 



For the proof of (i), from (3j, it follows that PS(i)'s 
maximize Ylk=i ^ (hoPs(k)) subject to the constraint d52l i: 
in other words, for the case of Pg(i) < P|(i), this value will 
be strictly decreased. Thus, (i) is proved. 

For the proof of (ii), it is equivalent to prove that 
ZLif(Pm) < EjL i < 3 < N. When 
j = i, by the assumption that Pg(i) < Pg(i), it follows 
*at ELi/(^W) < EUif(Pm)- For j = i + 
h if ELi/(^( fc )) > ELi/(^s( fc ))' it follows that 
f(PsU)) > f(P§U)) and thus P§(j) > P|(i) since f(x) is 
monotonically increasing. For Pg(i)'s, denote the index of the 
next source energy exhausted block after the j-th block as j s . It 
is then obtained that P|(fc) > P£(j) > pf(j) = P s {k), j < 
k < j s , since by Lemma ICTl Pg(i)'s are non-decreasing over 
i. Moreover, by O, since f(x) is concave, it can be shown 
that with the same energy budget over the first to the j-th 
blocks, if El =1 > Ei=i f(P*s(k)l " Allows that 

Efc=i p s( k ) > Ei=i p s( k )- Then, it is easy to check that 
Pg(j)'s violate the source energy constraint at the j s -th block, 
i.e., P£l=i P§(k) > PEi=i = EiLi Thus, 

it is obtained that ELi/(^sW) < Efc=i/(^sW) for 
j = i + 1. By using the mathematical induction method, 

it can be shown that ELi /(^s( fc )) < El=i f(P*s( k )^ 
1 < j < N, which suggests that (ii) is true. 

From the above analysis, it is proved that for the case 
of < ho < 1, P|(£) = Pj}{i)> i = I," - ,N, and thus 
P|(i)'s are optimal for Problem (P2). For the case of ho = 0, 
since there always exists one optimal source power profile 
that is non-decreasing over i (by Lemma ICTT i, without loss of 
generality, we can assume that Pg(i)'s are still non-decreasing 
over i. With the same argument as the case of < ho < 1, 
it is observed that claim (ii) is still true, which means that 
the maximum value with P|(i)'s is not larger than that with 
Pg(i)'s. In conclusion, this proposition is proved. 

Appendix D 
Proof of Proposition |53] 

Denote Pg(£)'s and P^(i + l)'s as the optimal solution for 
Problem (PI), Pg(i)'s and P£(i + l)'s as the optimal solution 
for Problem (P2). First, we prove that the condition given in 
Proposition 15.31 is sufficient. Assume that C(P^(k + 1)) < 
C{P s ( k )) ~ C{h P s {k)), 1 < k < i, and C(P^{i + 2)) > 
C(pf(i+l))-C(hoPs(i+l))- Due to constraint ([23), without 
loss of generality, we further assume that C(P^(i + 1)) < 
C(P<?W) -C(h Ps(i))- Then, for P£(i)'s and P*(i+1)% it 
is easy to check that only Scenario II or III given in Proposition 
14.21 can happen at the i-th block, which is further discussed as 
follows: 

1) Scenario III happens for the z-th source message in 
Problem (PI), i.e., P*{i) = P , P*(i + 1) = P + ^ - 1, 
and thus the source energy constraint at the i-th block is 
not satisfied. Then, there exists 5 > 0, such that the newly 
defined source power allocation Ps(i) = Pq + 6, Ps(i + 1) = 
Pc + — 1 — 5 satisfies the source energy constraint at the i- 
th block. Moreover, as for the NDC case, we can increase 
the binning rate of i-th source message with the amount 
C(P|(i + 1)) - C(P s (i + 1)) (note that this operation is 



15 



possible since this amount is less than C(Ps(i)) — C(PJ(i))), 
and transmit it at the (i + l)-th relay transmission. With the 
above scheme, it is easy to check that the sum rate over these 
two blocks is strictly improved. 

2) Scenario II happens at the i-th block for Problem (PI). 
Consider two subcases: (a) If the relay energy is not exhausted 
at the end of each TV-block transmission, as for the NDC case, 
we can increase Rs(i + 1), and use the available relay energy 
at the (7V+l)-th block to transmit the increased part in Rb{i+ 

1) , which strictly improves the throughput of the DC case; (b) 
If the relay energy is exhausted, there must exists k > i + 2, 
such that P R (k) > P^(i + 2). Thus, there exists < 5 < 
- p j?.( fc )-^ l H( i + 2 ) _ Define a new re i a y power allocation satisfying 
the relay energy constraint as Pr(i + 2) = P R (i + 2) + 5 and 
Pii(k) = Pfi(k) — S. By increasing the binning rate of the 
i-th message with the amount C(P R (i + 2)) - C(P R (i + 2)) 
and decreasing that of the fc-th message with C(P^(k)) — 
C(P R (k)), it can be shown that the new scheme improves the 
sum rate of these two blocks, since C(P R (i + 2)) — C(P^(i + 

2) ) > C{P^{k))-C{P R {k)). Based on 1) and 2), the "if" part 
is proved. 

Next, we prove that the condition in Proposition 15.31 is 
also necessary. Assume that C(P R (i + 1)) < C(PJ(£)) - 
C{h Pg{i)), Vi £ {1, • • • ,N}. It is easy to check that with 
the power allocation PJ(i)'s and P R (i + l)'s of Problem (P2), 
the achievable rate of the i-th block for Problem (PI) is given 
as R(i) = min{C(P|(i)),C(J^(* + l))+C(ftoP|(i))} = 
C (P R (i + 1)) + C (h Pg{i)), which equals that of the NDC 
case. By searching over the whole feasible set, the throughput 
of the DC case will be no smaller than that of the NDC 
case. Together with Proposition 13.21 it is obtained that the 
throughput of the two cases are identical, and thus the "only 
if" part is proved. 
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